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Abstract 

The dynamics of cosmic scalar fields with flat potential is studied. Their contribution to the 
expansion rate of the universe is analyzed, and their behaviour in a simple model of phase 
transitions is discussed. 
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The universe is observed to expand Q, and this expansion is suspected to have 
accelerated at a late stage of its evolution ^ [J. Furthermore, observation of 
the cosmic microwave background indicates the universe is spatially flat Q. 
The energy density of scalar fields has been recognized to contribute to the 
expansion of the universe j|, ||, and has been proposed to explain inflation 
||, as well as the presently observed expansion || [K], 11, 13 1; in ref. [14] an 



attempt has been made to confront the data with the predictions for a minimally 
coupled scalar field with an a priori unknown potential. An overview of models 
with many references has been given in |Iq[ . 

As suggested by the standard model of the electroweak interactions and its 
many proposed small-distance extensions with or without supersymmetry, the 
dynamics of the universe involves a large number of scalar fields interacting 
with gravity, gauge fields and fermionic matter. The observed isotropy and 
homogeneity of the universe do not allow for the existence of long-range electric 
and magnetic fields, but neutral scalar fields can have non-trivial dynamics in 
an expanding FRW-type universe. 

For a number of minimally coupled scalar fields (pi with a potential V[<j>\ the 
relevant dynamical equations for a flat universe (k = 0) are 

i 

where H is the Hubble parameter, and p is the contribution of matter to the 
density of the universe. It should be noted, that we absorb any cosmological 
constant in a constant contribution to the potential V [</>]. 

Constant scalar fields. As observed in j|, ^, [?J simple non-trivial solutions to 
the coupled equations ([!]) are provided by constant fields at stationary points 
of the potential: 

Vi = 0, V m + p=^. (2) 

Here V m is the value of the potential at the stationary point (presumably a 
minimum), and p represents the density of ordinary relativistic and cold non- 
relativistic matter, with present densities p r q and p nT Q. Then the last equation 
becomes in terms of the scale factor o~(t) = a(t)/ao: 

•2 ,2 2 , 87r C ( Prfl . Pnr,o\ ,2 8irGV m fo s 

a =ha +—1^+—), h=^-, (3) 

provided V m > 0. Hence for large universes the late-time behaviour of the scale 
factor is described by 

o(t) ~ a e h ^- t0 \ (4) 

This is consistent with the observed present state of the universe if the Hubble 
parameter is h ~ 70 km/sec/Mpc, of the order of the inverse life time of the 
universe. As is well-known, the corresponding value of V m is extremely small: 
V m = p c — 1CT 122 in Planck units. 

In contrast, for V m = 0, we obtain the standard radiation/matter dominated 
solutions with a(t) ~ t K , where « = 1/2 or k = 2/3. If V rn = -3w 2 /(8vrG) < 0, 
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there are oscillating solutions with maximal amplitude a = aao given by the 
solution of 

2_4 8ttG . _ . . . 

id a — {pnrfl O + Pr,0) = 0. (5) 

The last two cases can be consistent with a flat universe only if the density of 
matter and radiation equals or exceeds the critical density p c , which is definitely 
contradicted by the best current estimates of the densities of luminous and cold 
dark matter. 

Flat potentials. A different evolution of the Hubble parameter H{t) is obtained 
by allowing one or more scalar fields to depend on time. We consider the 
situation where the first field is dynamical: fa = (p(t), whilst all other fields are 
constant: fa = v% (i = 2,3, ...,N). Then V(tp) = V[fa = ip;fa = vi\. Thus we 
have 

Cp + Wy + V' = Q, \^ + v( {p ) +p = ^L. ( 6 ) 

In general any solution depends sensitively on the initial value (fi n it, which is 
difficult to control. However, this problem is absent if the potential is invariant 
under field translations <p = (p + e. Therefore it is of some interest to study flat 
potentials V(<p) = Vq. Such potentials are characteristic for Goldstone bosons, 
and they occur in many supersymmetric models. From the flatness of V one 
derives a conservation law for the field momentum per co-moving volume: 

(p + = =>■ a 3 (p = 7 = constant. (7) 

After multiplication by <r 6 the other equation then becomes 



The solutions can be characterized qualitatively as in the case of constant scalar 
fields: 

a. For Vq > 0, the potential provides the leading term, and for large a we find 
as before: 
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a(t) ~a Q e h{ t- t0 \ h = ^j 



b. For Vq = 0, the leading term is provided by the matter density p, and the 
evolution of the scale factor is again of the form a(t) ~ t K , where k = 1/2 for 
relativistic matter, and k = 2/3 for cold non-relativistic matter. 

c. Finally, for Vo < there are oscillating solutions, with the maximal scale 
a = aao a solution of 

^ +pa 6 + V a 6 = 0, (10) 
which is a direct generalization of eq.(||). 

Matter- dominated regime. The precise evolution of the scale factor and the 
scalar field can be solved analytically in the matter-dominated regime, in which 
pa 3 = constant. In the following we describe the solutions in some detail. 
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It is convenient to use the dimensionless scale variable a = a/ciQ = 1/(1+ z). 
Assuming Vq ^ 0, eq.(@) then becomes 



'da r 
~dt 



2AttGVc 



a 3 + 



+ 



po v 

2Vo) ' 2V 4K 2 



(11) 



where p(a) = pojo 3 . As this equation is invariant under time reversal, at any 
time we can have an increasing and a decreasing solution; we always exhibit 
here the solution with growing scale factor at earely times. For positive values 
of Vq we must distinguish three cases: 
a. 7 2 > Po/2V > 0. Then the solution of eq.(|ll|) is 



a(t) 



o(t) 
a 



A sinh3/i(t - t - A) 



Po 
2V 



nl/3 



with 



2V 4K 2 ' 



3hA 



-arcsinh 



2^0/ 



and /i given in eq.(^). The corresponding solution for the scalar field is 



1 



po 



arcsinh 



Vl2irG I 

with the intergration constant given by 



J_ 7 + Po ^ 

n-3 







(12) 
(13) 

(14) 



7 +Po 

2i 2 v -pl 



f2 = arcsinh 



As all equations are invariant under reflections of the scalar field: p 
there is another solution 



<p(v) = Po + 



-J_(. 



arcsinh 



(7 



1_ 7 2 + po cr 3 

3 



27 2 Vb - P 2 , 



n 



(15) 



-<P, 



(16) 



b. 7 2 = Pq/2Vo. In this case the constant terms in the square brackets in (|l 
cancel, with the result that 



a(t) 



1+1 + 



Po 
2V 



,3h(t-t Q ) 



1/3 



The solution for the scalar field reads 

1 



tp{a) 



Po + 



v / 12^G 



hi 



J_ 7 2 + po & 
a 3 7 2 + po 



(17) 



(18) 



To obtain the coefficients in front of and inside the logarithm, we have used the 
relation between Vq and h, and the special relation between Vo and po and 7. 
c. 7 2 < Pq/2Vq. The solution of eq.(ll) becomes 



a{t) 



A cosh[3/ l (t - t - A)] - i£_ 



1/3 



(19) 



3 



with 



A 



3hA = arccosh 



1 

A 



1 + 



Po 
2V 



The scalar field evolution is given by 



<p(a) = ifo =F 



1 



( 



arccosh 



_L_ 7 2 + po a 3 







(20) 



(21) 



where the constant is 



= arccosh 



7 + Po 
Pi ~ 2^o7 2 



(22) 



Observe, that in all three solutions above if diverges for a — > 0, as does 9?. 
However, long before this limit is reached the universe becomes relativistic, and 
our model must be modified. The initial conditions for the model discussed 
here should rather be set at the onset of the matter-dominated aera. 
d. If Vq = equation (@) implies 



a(t) 



1 + JUttG (2p + 7 2 ) (* - t ) + 6nGp (t - t y 



1/3 



(23) 



(We have chosen the solution which can vanish at t = 0.) This is the usual 
cold matter dominated universe. The scalar field equation is most easily solved 
directly in terms of time: 



ip(t) = ip Q =F 



arctanh 



s/VbvG p (t - to + A) 
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with 



A 



/ 7 2 + 2p 



This solution is asymptotically constant; 



Q, = arctanh W 1 + 



2/50 



</>0 ± 



1 



( 



arctanh 



'1 + 



2po 

7 2 



(24) 



(25) 



(26) 



The field can have reached its asymptotic value at present only if <C po- 
e. For Vq < 0, the scale factor takes the periodic form 



with 



a(t) 



A 



Po 
2\V \ 



+ A sin3u(t - t - A) 



1/3 



Po + 7 2 



4V 2 2|Vb|' 



'8vrG|y | 



(27) 



(28) 
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and 

3ujA = arcsin 
The scalar field solution reads 

1 



Po 

A \2\V \ 



ip(a) = <p =F 



( 



arccosh 



cr 



7 2 + ppcr 3 



n 



where 



Q = arccosh 



T + Po 



(29) 



(30) 



(31) 



Actually, this is formally the same solution as for the scalar field in case (c). 
However, it evolves in a universe with a maximum expansion and a finite life 
time, rather than an infinitely expanding one. 

Equation of state. As cold matter is pressureless, the total pressure of the 
cosmic fluid comes from the scalar field: 



Ptot = Ps 



The total density is 



7 

2a 6 



7 



-V . 



Plot =Ps+P=7^ + ^ + V . 



It follows that the state parameter is given by 

Ptot 2V a 6 



Wtot 



7 



Ptot 



2Vo<7 + 2p a 3 + 7 2 



(32) 



(33) 



(34) 



For infinitely expanding universes with V > 0, in the large-scale limit the 
parameter becomes Woo = — 1- Inflationary behaviour sets in when wtot < —1/3. 
This happens at a critical scale a e given by 



0". 



Po 
4F 



1 + W1 + 



167 2 Vb 



Pi 



(35) 



Depending on the relative values of Vq and Pq/2~/ 2 , we have 

7 2 > P 2 o/2V 

7 2 = P 2 o/2V 

7 2 < P 2 o/2V 



cr 3 > 

a e > Vb' 



3 = Po. 

e Vb' 



(36) 



For vanishing potential Vq = the parameter is always positive, but approaches 
wtot i in the large-scale limit. For Vq < the state parameter wtot is always 
positive and becomes infinite at the turning point a = a. 
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The ratio of dark energy to cold matter is given by 

Ps 7 2 1 V 3 

y = — = „ ^ + — ° ■ ( 37 ) 



For Vo > it never becomes small, but it reaches a minimum for = 7 2 /2Vo, 
when 

y m = \ — — • (38) 
V Po 

Thus a situation in which y m < 1 can be reached only if 7 2 < Pq/2Vo. In this 
case 

< 7- ^ 

From this analysis we observe, that a recent acceleration of the universe con- 
jectured on the basis of the type la supernova data can be incorporated in 
these simple models for Vo > 0. By construction the development of the scalar 
field does not depend strongly on the initial conditions, as in more sophisticated 
models with tracking behaviour [16]. Also, we can arrange for the scalar density 



to be comparatively low in the early universe, and have accelerated expansion 
at a later time if we take the parameters in the domain < 7 2 < Pq/2Vq. 

Piece-wise constant potentials. The analysis of flat potentials can be extended 
in a straightforward way to the case of piece- wise constant potentials. The 
main new feature is that whenever the field reaches a potential step, we must 
determine and match its behaviour before and after it hits the discontinuity. 
This is done by rewriting the field equation in the form 

5(^) = ^ - (40) 

The second equation can be integrated under two different conditions: 

(i) if the field crosses the potential step (transmission), we obtain 

A7 2 

2^ = -AV, (41) 

with AV = V2 — V1 the magnitude of the potential step, a the scale at which the 
field reaches the discontinuity in the potential, and A7 2 = 7| — 7i the change 
in the square of the field momentum; 

(ii) if the field does not cross the potential step (reflection), the value of the 
potential does not change, and 7 2 is equal before and after the reflection. As 
the field now must run in the other direction as a function of time, we deduce 
that 72 = -71. 

The condition (^l|) can also be written in the form 

£ S + V 1 = § S + V 2 . (42) 

As matter is not created nor destroyed, this shows, that H is continuous during 
the change in the scalar field. In contrast, during transmission across the po- 
tential step the derivative of the field <p, the derivative of the Hubble parameter 
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H and the potential all change by a finite amount. This can also be derived by 
taking a derivative of the second eq. @ , which leads to the equations 

h i 2 1 m 2 + H Tr i 



8nG 2 r 2 r ' 8vrG 2 

From the continuity of p and H, and the finite discontinuity of V, the finite 
discontinuous changes in H and 93 are obtained directly: 

The same equations actually also hold trivially at reflection, as all changes in 
these quantities then vanish. In this case continuity of the field is obtained by 
taking a solution y?(<r) = <Pi(<j) before the reflection (a < a), and <p(cr) = ^2(0") 
after the reflection {a > a), with the two solutions related by 

ip 2 (a) = 2if - ipi((r), (p = ipi(a) = (p 2 (a). (45) 

Observe, that if the initial condition is such that 

^ + V 1 <V 2 , (46) 

then only the reflected solution exists. 

Obviously, the piece-wise constant potentials discussed here can be used 
as approximations for a variety of smoothly changing potentials. It is to be 
noted that the translation symmetry in field space is then lost, but one still has 
relative freedom in choosing initial conditions in any flat region of the potential. 

Phase transitions. In the course of evolution of the universe, phase transitions 
accompanied by symmetry breaking can change the value of the potential. We 
consider the simple situation that the potential changes from a value V\ before 
the transition time t c to a value V 2 after t c . For simplicity we take the transition 
to occur instantaneously. We can then treat the phase transition as a potential 
step occurring at the time t c . Eq.(|42"|) then implies, that during the transition 

where a c = o~(t c ). Note, that in this case we only have the possibility of 
transmission, not reflection of the field. As the matter density and Hubble 
parameter are not affected during the transition, the continuity of a{t) and 
(p(t) allow us to express the parameters of the solution for t > t c in terms of 
those of the solution for t <t r . 



In this context, consider again the parameter y m = y 27 2 V/ the minimal 
ratio of scalar to matter density. From eq.(^7|) it follows, that after a phase 
transition 

FO z Po a c 
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Thus, whether the transition leads to a situation with larger or smaller y m 
depends both on the difference in potential energy density as well as on the 
difference in kinetic energy density of the field. In particular, the value remains 
the same if 

Vx = ll JL = V 2 ^ = V 1 (49) 

i.e. if the values of kinetic and potential energy are interchanged during the 
transition. In all cases where this condition is not met, the value of y m changes. 
If the potential energy density dominates over the kinetic energy density, than 
lowering V will in general lower y m , whereas in kinetic-energy dominated situ- 
ations a lowering of V will give rise to a larger value of y m after the transition. 
This follows from the observation, that eq.(|48|) can be replaced by the equivalent 
expression 

V 2 m2 -y 2 ml = 2 jt M " V2) U + V 2 - . (50) 

As generically the nature of the solution for a(t) and ip(t) depends on whether 
y^ is positive, zero or negative, and in the first case on whether it is equal 
to, larger or less than unity, the solutions before and after a phase transition 
may belong to different classes a-e discussed above. It should be noted that 
the analysis also applies to the case where an originally constant field (71 = 0) 
becomes dynamical (72 / 0), in which case V\ = V 2 + 7!/2of . 

In conclusion, massless scalar fields associated with flat potentials display 
simple but interesting cosmological dynamics. As they are associated with 
Goldstone bosons, they can play a role in phase transitions during which they 
change their dynamics. Finally, it has been shown how to piece together the 
solutions for potentials with discrete jumps, which under certain conditions can 
give qualitative insight in the behaviour of non-constant potentials. 
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